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Abstract. In this paper we study the large-time behavior of classical solu- 
tions to the two-species Vlaaov-Maxwell-Boltzmann system in the whole space 
. The existence of global in time nearby Maxwellian solutions is known from 
1371 in 2006. However the asymptotic behavior of these solutions has been a 
challenging open problem. Building on our previous work 1121 on time decay for 
the simpler Vlasov-Poisson-Boltzmann system, we prove that these solutions 
converge to the global Maxwellian with the optimal decay rate of 0(t~ 2 + 2F) 
in L|(LJ)-norm for any 2 < r < 00 if initial perturbation is smooth enough 
and decays in space-velocity fast enough at infinity. Moreover, some explicit 
rates for the electromagnetic field tending to zero are also provided. 
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1. Introduction 

The Vlasov-MaxweU-Bohzmann system is an important physical model to de- 
scribe the time evolution of dilute charged particles (e.g. electrons and ions) under 
the influence of their self-consistent internally generated Lorentz forces [33] • The 
existence of global in time solutions to this Cauchy problem which are perturba- 
tions of a Maxwellian equilibria is known since 2006 in 37 . The time rate of 
convergence to equilibrium is an important topic in the mathematical theory of the 
physical world. It has however remained an open problem to determine the large 
time behavior of these solutions. In this work, we resolve this problem proving con- 
vergence to Maxwellian with the optimal rate of 0(t~2+^7) in L|(_L^)-norm with 
2 < r < c» in three spatial dimensions. 
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1.1. Equation and reformulation. For two species of particles, the Vlasov- 
Maxwell-Boltzmann system is given in its non-dimensional version as 

dtf++^- V,/+ + {E + ^xB)- Ve/+ = 0(/+, /+) + Q{f+J^), 

dtf- + e • V,/_ -{E + ^xB)- V^/_ = Q{f-J+) + QU-J-). 

It is coupled with the Maxwell system 

(2) dtE-V,xB = - f af+-f-)d^, 

(3) dtB + V,xE = 0, 

(4) V.-E^ f {f+~f^)d^, V, B = 0. 
The initial data in this system is given as 

(5) /±(0,x,O =/o,±(a;,0, E{0,x) ^ Eo{x), B{0,x) ^ Bo{x). 
The initial data should satisfy the compatibility conditions 

y^-Eo^ [ (/o,+ -/o.-)de, V,-Bo = 0. 

Here the unknowns are f± ~ f±{t,x,£,) > and E{t,x), B{t,x), standing for 
the number densities of ions (+) and electrons (— ) which have position x = 
{xi,X2,X3) e M.^ and velocity ~ {^i,£,2,£.3) & at time t > 0, and the elec- 
tromagnetic field, respectively. Q is the bilinear Boltzmann collision operator for 
the hard-sphere model defined by 

f^f{t,x,0, f' = fit,x,i'), =.9(t,x,e*), 9:=9it,x,C), 

We will study solutions to this system which are initially perturbations of the 
Maxwellian equilibrium states. 

We write the normalized global Maxwellian as 

M = M(e) = (27r)-3/2e-l«l'/2. 

Then we linearize around the perturbation u in the standard way 

Denote the column vector [•,•] as follows: / = [/+,/-], u = Then the 

Cauchy problem ©, 0, ©, g]), (0 can be reformulated as 

(6) dtu + ^- V^u + q(E + ^xB)- V^u - E ■ ^M^^^qi 

= Lu + ^E ■ -I- r(u, u), 

(7) dtE^V^x B = - I £m^/^{u+ - u^)dt 

(8) dtB + \7^xE = Q, 

(9) V^-E^ j M.^''^{u+-u^)di, V^-B = 0, 
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with initial data 

(10) u±iO,x,0=uo,±ix,0, EiO,x)=Eoix), B{%x) = Bo{x), 
satisfying the compatibihty condition 

V^-Eq^ [ M'/^{uo.+ -uo,-)d^, V,-Bo=0. 

JS.3 

Here, q = diag(l,— 1), qi = [1,-1], and the linearized collision term Lu and the 
nonlinear collision T(u,u) are respectively defined by 

Lu — [L+u, L_u], r(u, v) = [r+(u, w), r_(M, v)], 

with 

In the sequel we will study the asymptotic behavior of solutions to this system. 

For the linearized collision operator L, one has the following standard facts [3]. 
L can be split as Lm — —iy(^)u + Ku, where the collision frequency is given by 

Notice that v{£) ^ (1 + The null space of L is given by 

U = span {[1, 0]Ml/^ [0, l]M''\ K„ ^^]M^'^ (1 < z < 3), m\ \e]m^'^] . 

The linearized collision operator L is non-positive and further — L is known to be 
locally coercive in the sense that there is a constant Aq > such that [3]: 

(11) -/ uLud^>Xof i^mil-Fjul^dt 

JR3 JR3 

where, for fixed {t, x), P denotes the orthogonal projection from L| to J\f. 

1.2. Notations and main results. Throughout this paper, C denotes some pos- 
itive (generally large) constant and A denotes some positive (generally small) con- 
stant, where both C and A may take different values in different places. In addition, 
A ^ B means XA < i? < -^A for a generic constant < A < 1. For any integer 
m > 0, we use i/", iJ™, to denote the usual Hilbert spaces i?™(R; x R^), 
-ff™(R"), iJ™(]R^), respectively, where L^, L^, L| are used for the case when m = 0. 
For a Banach space X, \\ ■ \\x denotes the corresponding norm, while || • || always 
denotes the norm || • ||^2 or || • 11^2 for simplicity. We use (•, •) to denote the inner 
product over the Hilbert space L|, i.e. 



and for simplicity, (•, •) is also used as the inner product over when there is no 
possibility of confusion. For r > 1, we also define the standard time-space mixed 
Lebesgue space Zr = i|(i^) = ^^(R^; i''(Mi^)) with the norm 

2/r \ 1/2 

Iffb.-I/ (/ \9i^,0rdx\ dA , g = gix,0(^Zr. 
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For multi- indices a — [ao, ai, ck2, as] and /3 — [/3i , /32 , /^a] , we denote 

QQ _ aaoaai aa2 «Q!3 a/3i a/32 afe 
C/3 - C* Caja "^3 C*^! ^'^^ ■ 

The length of a is |a| = ao + ai + q;2 + aa and the length of /? is |/3| = /3i + /32 + /Ss- 
For simplicity, we also use dj to denote d^^ for each j = 1,2,3. 

For an integrable function (7 : M'' — > M, its Fourier transform is defined by 

f 3 
g{k) = Tg{k) = / e-^''"-''g{x)dx, x • fc = V xjkj, k £ R^ 

where z = \/— T G C is the imaginary unit. For two complex vectors a,b G C^, 
(a I b) — a ■ b denotes the dot product over the complex field, where b is the 
ordinary complex conjugate of b. Note further that we will use a subscript j in an 
equation number to denote the j-th equation displayed. Precisely, for example (fl])9 
refers to the second equation in ([T|). 

Given a solution [u{t,x,£^),E{t,x),B{t,x)] to the Vlasov-Maxwell-Boltzmann 
system ((6|)- (|10p . an instant full energy functional is defined as a continuous function 
£N,m{t) which satisfies (IT^ . The instant high-order energy functional £^}f „^{t) is 
defined similarly, satisfying instead (1131) . In particular we use 



(12) sna^)-^ E E wim.Bimw 

\a\ + \P\<N \a\<N 

^N,n.{t)- E lk^5"^t)||2+ ^ ||^f9|{I-PHi)f 

l<|Q|<Af \a\ + \fi\<N 

(13) + E W[E{t),B{t)]r + \\E{t)\\\ 

l<\a\<N 

and we define the dissipation rate I'Ar,m(i) as 

l<|Q|<Ar |q| + |;3|<7V 

(14) + E \mE{t),B{t)]r + \\E{t)\\\ 

1<|q|<JV-1 

for integers N and m. For brevity, we write fAr(t) = £N,o{t), ^Ar(0 = '^]v.o(0 and 
'DNit) = T^N,o{i) when m = 0. Throughout this paper we assume TV > 4. 

From [37] and a little additional efforts about estimates of the electromagnetic 
field [E{t, x), B{t, x)], we have the following existence result. 

Proposition 1.1. Let [uq, Eq, Bq] satisfy © and 

/o,± (x, - M + Ml/ V±(2^, > 0. 

There exists £N{i) such that if £n{0) is sufficiently small, then the Cauchy prob- 
lem ([7]), ([5]), (O, and pn]) to the Vlasov-Maxwell-Boltzmann system admits a 
unique global solution [u{t,x,$^), E{t,x), B{t,x)] with 

/± (t, X, = M + Ml/ V (i, X, > 0, 
and the Lyapunov inequality 

(15) ^fjv(i) + A2?Ar(0 <0, 
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for any t > 0. 

Moreover, on the basis of the above existence resuh, one can obtain some 
velocity-weighted or high-order energy inequahties stated as fohows. 

Theorem 1.2. Suppose that all the assumptions of Proposition [177] hold. For any 
given m — 0,1,2,..., there exists £N,m{t) such that «/ fjv,m-i(0) is sufficiently 
small and f iv.m(O) is finite, then 

(16) ^fjV.mW+AI^A'.mW <0, 

at 

holds for any t > 0, where fjv,-i(0) = £n{Q) is set and A may depend on m. In 
addition, if £^(0) is sufficiently small, then there is f]v(0 = oi^) such that 

(17) j^SUt) + APiv(t) < Ci\\iy'^^V,Pu{t)f + \\V, x B{t)f), 
holds for any t > 0. 

From the proof later on, the above nonlinear energy estimates together with 
the time-decay estimates on the linearized non-homogeneous system indeed lead to 
optimal time-decay rates of the instant energy functionals EN.mif) m = 0, 1, 2, • • • 
and S^{t) under some additional regularity and integrability conditions on initial 
data. Precisely, our main results in this paper are stated as follows. Set „i as 

(18) ej,™ = fj,m(0) + ||«olll, + ||[i?o,So]|lii, 
for integers j and m. 

Theorem 1.3. Let [u{t, x, ^), E{t, x), B{t, x)] be the solution to the Cauchy problem 
(IS]), ([HI), and pO|) of the Vlasov-Maxwell-Boltzmann system obtained in 
Proposition For any fixed m = 0,1,2,..., if ^N+2,mvi is sufficiently small 

where m V 1 == max{m, 1}, then 

(19) £NMt)<CeN+2,rn{l+t)-i, 

holds for any t >0. In addition, if eN+5.i is sufficiently small then 

(20) £Nit)<C^N+5Al + trK 

holds for any t > 0. 

Furthermore, we have the optimal decay rate in the L'" norm for any 2 < r < cx3: 

Corollary 1.4. Let [u{t, x, ^), E(t, x), B(t, x)] be the solution to the Cauchy problem 
0, and (|10|) of the Vlasov-Maxwell-Boltzmann system obtained in 

Proposition Suppose that eg^i and eg, 2 are sufficiently small, then for any 

2 < r < 00, we have the following estimates for any t > 0: 

(21) \\u{t)\\z^ + \\B{t)\\L'-<C{l + t)-i+^, 

(22) ||{I - P}u{t)\\z,. + ||([1, -l]Mi/2,«(t))|U. + \\E{t)U. < C{1 -f tri+^. 

Remark 1.5. From Theorem l3.11 which proves the time decay rates of the linearized 
system, the solution [u, E, B] to the linearized homogeneous Vlasov-Maxwell-Boltzmann 
system satisfies 

(23) |ivx<)|| + ||v™[i?(<),i?(t)]|| <c(i + i)-i-^. 
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for any t > 0, where to > is an integer, C depends on the initial data, and the 
initial data is smooth enough and decays fast enough in x, ^ at infinity. By applying 
further the optimal Sobolev inequality as in (|122p and interpolation, implies 

\\u{t)\\z,. + \\m),Bit)]Ur^<Cil + tri+*= 

for any t > 0, where 2 < r < oo. Thus, for the nonlinear system, the time decay 
rate in (j21l) for u and B is the same as the above one in the linearized case, and 
the faster time decay rate in for {I - P}u, ([1, -IJM^/^, u{t)) and E{t) is due 
to the special structure of the Vlasov-Maxwell-Boltzmann system. On the other 
hand, when the electromagnetic field disappears, the spectral analysis as in |46j for 
the one-species Boltzmann equation in the hard sphere case shows 

l|Vre*(-«-^^+^)uo|| < C(l + i)-i(t-^)-^(ho||z, + ||V>o|l) 

for any t > 0, where to > is an integer and 1 < p < 2. Thus, the time decay rate 
of u in (PT|) given by 

Mt)\\z^ <C{l + t)-i+^, 2<r<oo, 

is optimal in the sense that it is the same as in the cases of the Boltzmann equation 
without any force and of the linearized Vlasov-Maxwell-Boltzmann system. 

Notice here that the electric field decays faster than the magnetic field in with 
2 < r < oo. Indeed recall that the pure homogeneous Maxwell's equations usually 
conserves the energy. Thus the decay of these terms is truly a non-linear effect 
which results from the coupling with the Boltzmann equation. The same feature 
has been observed in [8] for a study of the Euler-Maxwell system with relaxation, 
where the Green's function of the linearized system is presented in detail and thus 
leads to the optimal large time behavior of each component in the solution. 

1.3. Literature, new ideas, and future possibilities. It was pointed out by 
Villani in [?7] , that there exist general structures in which the interaction between 
a conservative part and a degenerate dissipative part will lead to convergence to 
equilibrium; this property has been called hypocoercivity. This paper provides a 
concrete example of hypocoercivity for the nonlinear two-species Vlasov-Maxwell- 
Boltzmann system in the framework of perturbations. We notice that the general 
theoretical framework in [IT] can not be directly applied here. On the other hand, 
we hope to further develop the hypocoercivity theory in the future to include some 
degenerately dissipative kinetic equation coupled with a class of hyperbolic systems. 

There has been extensive investigations on the rate of convergence for the non- 
linear Boltzmann equation or related spatially non-homogeneous kinetic equations 
with relaxations. In what follows let us mention some of them. In the context of 
perturbed solutions, the first result was given by Ukai [45], where the spectral gap 
analysis was used to obtain the exponential rates for the Boltzmann equation with 
hard potentials on torus. The results in [45j were improved by Ukai- Yang |46j in 
order to consider existence of time-periodic states in the presence of time-periodic 
sources, which was later extended by Duan-Ukai- Yang-Zhao [TC] to the case with 
time-periodic external forcing by using the energy-spectrum method; see also |15j 
for a summary of general applications of the energy-spectrum method when some 
given small external forcing occurs. Recently, Strain-Guo |42j developed a weighted 
energy method to get the exponential rate of convergence for the Boltzmann equa- 
tion and Landau equation with soft potentials on the torus. Earlier but along the 
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same line of research, Strain-Guo [41j developed a general theory of polynomial 
decay rates up to any order in a unified framework and applied it to four kinetic 
equations, the Vlasov-Maxwell-Boltzmann System, the relativistic Landau-Maxwell 
System, the Boltzmann equation with cutoff soft-potentials and the Landau equa- 
tion all on the torus. 

Another tool is entropy method which has general applications in the exis- 
tence theory for nonlinear equations. By using this method as well as the elabo- 
rate analysis of functional inequalities, time-derivative estimates and interpolation, 
Desvillettes-Villani [4 obtained first the almost exponential rate of convergence 
of solutions to the Boltzmann equation on torus with soft potentials for large ini- 
tial data under the additional regularity conditions that all the moments of / are 
uniformly bounded in time and / is bounded in all Sobolev spaces uniformly in 
time. See Villani [4^ for extension and simplification of results in [4 still con- 
ditionally to smoothness bounds by further designing a new auxiliary functional. 
Notice that |41| provided a very simple proof of [4] for the unconditional pertur- 
bative regime. Recently, by finding some proper Lyapunov functional defined over 
the Hilbert space, Mouhot-Neumann [35] obtained the exponential rates of con- 
vergence for some kinetic models with general structures in the case of torus; see 
also [17] for the general study. An extension of [3S] to models with a single mass 
conservation law in the presence of some confining potential forces was given by 
Dolbeault-Mouhot-Schmeiser |7j. Recently, Duan [TTj developed a general method 
to deal with a class of linear degenerately dissipative kinetic equations with or 
without confining forces even when several physical conservation laws are present, 
specifically the linear Boltzmann equation with a parabolic confining force included. 
We remark that some ideas in for finding the dissipation of all the macroscopic 
components of the kinetic equation will be essentially used in the proof of our main 
results of this paper. 

Besides those methods mentioned above for the study of rates of convergence, 
the method of Green's functions was also founded by Liu-Yu [32J to expose the 
pointwise large-time behavior of solutions to the Boltzmann equation in the full 
space M^. 

In addition, concerning the Vlasov-Poisson-Boltzmann system, Glassey-Strauss 
[l8] studied the essential spectra of the solution operator and they [17^ also used 
the method of thirteen moments by Kawashima 29 to obtain the exponential time- 
decay rates of solutions on torus. A slow time-decay rate for the solution in the 
whole space was obtained by Yang-Zhao '49| on the basis of the pure energy method 
and a time differential inequality. Notice the series of works [9] [TOl [12 investigating 
solution spaces without initial layer. Then Duan- Strain [12] obtained the optimal 
time-decay rates of solutions to the one-species Vlasov-Poisson-Boltzmann system 
in M'^. The key of the method to study hypocoercivity provided by [12] is to carefully 
capture the full dissipation of the perturbed macroscopic system of equations with 
the hyperbolic-parabolic structure, which is in the same spirit of the Kawashima's 
work [28 . For the linearized time-decay analysis, instead of using the compensa- 
tion function as in [28) [29] , the main idea of [12l [H] is to design some interactive 
functionals in order to take care of the dissipation of the degenerate part in the 
solution. 

Here, we mention that if there is no coUisional effects as in the Vlasov-Poisson 
system, then the so-called Landau damping comes out. This was recently studied by 
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Mouhot-Villani [36 on torus, where it was shown that even though in the absence 
of kinetic relaxation, in the analytic regime the solution still converges weakly to 
large-time states determined by the initial data and the nonlinear system itself. 

For the more intricate Vlasov-Maxwell-Boltzmann system, which includes the 
hyperbolic coupling, as mentioned before, rapid polynomial time decay up to any 
order on the torus (for the solutions from [24) ) was shown in Strain-Guo [41) . 
Different from the case of torus, solutions in the whole space are dispersive and 
hence only the slower polynomial time decay is expected as in [T^]. Moreover, 
because of the weaker dissipative property of the Maxwell equations than in the 
case of the Poisson equation [T2j, it is difficult to study the large-time behavior of 
solutions to the Vlasov-Maxwell-Boltzmann system in the whole space. 

Finally, we also mention some of the results on the existence theory of the simpler 
Vlasov-Poisson-Boltzmann system and related kinetic equations: global existence 
of renormalized weak solutions with large initial data was studied in (BJ [SI 134] . 
global existence of classical solutions near Maxwellians can be found in |31] [30l [48] , 
[23[ [21] [37] and (Qj [13] , and global existence of solutions near vacuum was shown in 
[22[ [14] . We anticipate that the time derivatives may be removed from the existence 
theory [37] and our main results herein using methods from [27[ [9] . 

To obtain the results from this paper in Theorem ll.2[ Theorem ll.3[ and Corollary 
11.41 the starting point is our first paper [12] and the existence theorem from [37] . 
There are additionally several new ingredients which we expect can be useful in 
a number of other contexts, some of which we list after explaining the key new 
elements in our proof just now. First of all, from the hyperbolic nature of Maxwell 
equations it seems to be necessary to lose a derivative in estimates of the dissipation 
rate for the electro-magnetic field. This makes it very difficult to prove time decay 
without paying the price of extra derivatives. This problem is also present in the 
case of the torus [41] , but it is much more problematic in the whole space because 
the dissipation is necessarily dramatically weaker [37 . To surmount this difficulty 
in the whole space, we first spend some time developing a new pointwise time- 
frequency Lyapunov functional which measures directly the weak dissipation (and 
spatial derivative loss) of the full Vlasov-Maxwell-Boltzmann system; this is given 
in Theorem 13.51 Here new estimates are developed to measure the macroscopic 
dissipation and the very weak electro-magnetic dissipation in the whole space by 
introducing two interactive time-frequency functionals, 

which are defined by ([55]) and ([57|. respectively. After this we can make use of 
the usual frequency splitting method in Section [3?3] in order to take advantage of 
the weak dissipation and prove the linear decay which is stated in Theorem 13.11 
We furthermore discuss a new time-frequency (time-derivative) splitting method 
in Remark 13.61 which is related to the general derivative loss phenomena. For the 
low frequency part, the solution decays like the heat kernel, while for the high 
frequency part, it decays in time with the algebraic rates of any order as long 
as initial data is regular enough in space variable. We are optimistic that this 
approach can be useful in other hyperbolic or mixed hyperbolic-parabolic systems 
where derivative loss is present. For the study of some macroscopic systems with 
this kind of regularity-loss property, we here mention the recent work [551 HSl HI] 
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by Kawashima and his collaborators. Another new point is to introduce the time- 
weighted method to handle the time-decay rates of the full instant energy functional, 
which seems necessary because of the regularity- loss property of the electromagnetic 
field. Furthermore, to obtain the optimal time decay rates in L"^ with 2 < r < oo 
for the separate components of the solution as stated in Corollary 11.41 we use an 
optimized Sobolev inequality as in (|122|) . 

To finish we mention some expected future applications of these methods. As 
a result of scaling and the Lorentz invariant Maxwell equations, a very physically 
relevant collision kernel to apply when the kinetic equation is coupled with its inter- 
nally generated electromagnetic forces is the relativistic Landau-Maxwell system, 
for which the existence theory is known from [40]. We expect that with a suitable 
relativistic modification of the moment equations in Section [21 then the methods 
in this paper and p2] can be used to observe that the all of results herein with 
the same conclusions can be obtained for the relativistic Landau-Maxwell system 
in the whole space. The details of this approach would be quite complicated and 
interesting to carry out. 

Additionally we expect that the methods developed in this paper and [12] can be 
useful in several other physical contexts. In particular one should be able to com- 
bine the methods from 41 with our methods here to obtain optimal convergence 
rates for the relativistic Boltzmann equation with soft potentials |38j in the whole 
space. More recently, a global existence theory has been developed for perturba- 
tive solutions to the Boltzmann equation without the angular cut-off assumption 
for the full range of inverse power-law potentials on the torus [SQJ [19] . The meth- 
ods used there yield exponential decay for hard and moderately soft potentials 
(7 -I- 2s > 0), as well as "almost exponential decay" for the softest potentials. The 
methods developed in this paper and [41j [12] have also been quite useful in obtain- 
ing the optimal decay rates in the whole space for the these solutions flQl [19] for 
all of the hard and soft-potentials [39] . Finally, we would like to mention that in 
the direction of the non cut-off Boltzmann equation, a study about the qualitative 
properties of classical solutions, precisely, the full regularization in all variables, 
uniqueness, non-negativity and the convergence rates to the equilibrium, has been 
also provided by Alexandre, Morimoto, Ukai, Xu, and Yang in [1] and [2] (and the 
references therein), where the optimal time decay rates for the hard potentials are 
obtained by using the compensation function method developed in 29 . 

1.4. Organization of the paper. In Section [21 we define the macroscopic pro- 
jector and derive some moment equations up to third-order. In Section [S] we prove 
Theorem 13.11 about the time-decay property of solutions to the linearized Vlasov- 
Maxwell-Boltzmann system with microscopic source terms. This is based on The- 
orem 13.51 for the construction of a time-frequency Lyapunov functional resulting 
from estimates on the microscopic dissipation, macroscopic dissipation and electro- 
magnetic dissipation proved in Section [3. 2 1 Then, the proof of Theorem [ST] is given 
in Section 13.31 In Section [4] we study some nonlinear energy estimates including 
the time evolution of the equivalent velocity- weighted instant total energy and the 
high-order energy to finish the proof of Proposition 11.11 and Theorem 11.21 Then in 
Section [5] we use the decay of the linearized solution with nonhomogeneous source 
terms together with those nonlinear energy estimates to bootstrap and obtain the 
nonlinear decay. In Section 15.11 and Section 15. 2| we prove Theorem 11.31 Finally in 
Section lO] we prove CoroUarv 11.41 
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2. Moment equations 



In this section, we begin with the representation of the macroscopic projector 
P and then derive from the perturbed system some macroscopic balance laws and 
high-order moment equations which are systems of first order hyperbolic equations 
of the macroscopic coefficient functions coupled with the high-order moment func- 
tions. 

Given any u{t,x,£,), one can write P in as 
(24) Pu = a+{t, x)[l, 0]Mi/2 + a^{t, x)[0, 1]M^/^ 

3 

since P is a projection from L| x to A/", where the coefficient functions a±(t, a:), 
b(t,x) = [bi{t,x),b2(t,x),b3{t,x)] and c(i,x) depend on u{t,x,£,). The expression 
((24| can then be rewritten as Pu — [P+u, P_it] with 



P±u = {a±{t, x) + b{t,x) ■ e + c{t, xM\^ - 3)}Mi/2_ 
Since the projection P is orthogonal we have 

which together with the form (IM)) of P imply 
a±=.(Mi/2,u±) = (Mi/2,P±u), 



(e.Mi/2,p±u), 



12 
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In the rest of this section we will derive the equations for these macroscopic variables 
and also the high-order moments as follows. 

First consider the following linearized system with a non-homogeneous source 
9 = [9+{t,x,0,9-{t,x,^)]: 



(25) 



f dtu± + e • V^u± T E ■ ^ L±M + g±, 

dtE~v.,xB = - I sm^/^{u+ - u^)dt 

dtB + V^xE^ 0, 

V^-E^l Mi/2(u+ - u_)df, W^-B = 0. 



Taking velocity integrations of (j25p i with respect to the velocity moments 

Ml/2, ^^Mi/2,j^ 1,2,3, i(|Cp-3)Mi/2, 

6 
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one has 



(26) dta± + V, • + V, • {I± - P±}u) = {M^^^, g±), 
dt[h + {I± - P±}«)] + a,(a± + 2c) T 

(27) +V, • (^e.Mi/2,{I± - P±}u) = (6Mi/2,g±+L±u), 

c+i((|ep-3)Mi/2,{I±-P±}t 



5, 



6 



(28) +iv, • ((|^|2 - 3)^Mi/2, {I± - P±}w) = 1((|C|2 - 3)M'/^g±+L±u), 

where we have set I = with I±u = u±. Define the high-order moment 

functions 0{u±) = (9i-,(u±))3x3 and A{u±) = (Ai(u±), A2(u±), A3(u±)) by 

(29) e,,(H±) = ((e,e,-l)Mi/2,7.±), A,(u±) = i^((|ep-5)e,Mi/2,w±). 

Further taking velocity integrations of with respect to the above high-order 

moments one has 



(30) 


^t[Q^^{{I± 




+ 2c\ + 2d^b,^eu{l±+g±), 




^te^Ji{I± 




+ d,h + d,b, + V, • (^Mi/^ {I± - P±}u) 


(31) 






-e,,(/±+.g±) + (Mi/2,5±), 


(32) 


dtM{l± - 


P±}U) ^ 


hSjC = Aj(?± +5±), 


where 








(33) 






-C- V:,{I±-P±}u + L±u. 



Here we used the moment values of the normalized global Maxwellian M: 

(1,M) = 1, (|e,f,M) = l, (|eP,M)=3, 
(|e,f|e™P,M)-l, j^m, 
(|e,|^M)=3, (|ene,f,M)=5. 

Additionally to derive pil) we have used 

In particular, for the nonlinear system (jH])-®; the non- homogeneous source g 
[g+{t, X, ^), g_ {t, X, ^)] takes the form of 

(34) g± = ±^E -^u^TiE + ^x B)- V^u± + T±iu, u). 

Then, it is straightforward to compute from integration by parts that 

(Mi/2,g±) =0, 
(^M^/^ g±) = ±Ea± ±bx B±{^M^/^, {I± - P±}u) x B 

hue - m'^^9±) = ±\b-E± \m^/^ {i± - p±}u) . e 

6 3 3 

+(i(|ep-3)Mi/2,r±(u,w)). 

6 
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Thus, the balance laws ([26l) - ([28| for the general case can be re- written as 



(35) 



(36) 



(37) 



dta± + V, • 6 + V, • i^M'/^, {I± - P±}u) = 0, 
dt[h + (6:Mi/2, _ p^i^^] ^ Q^^^^ ^ 2c) T 

+V, • (ee«Mi/2, {I± - P±}u) = ±E,a± ± [6 X B], 

±[(eMi/2,{I± -P±H X B], + (e.Mi/2,L±u + r±(u,u)), 



c+i((|Cp-3)Mi/2,{I±-P±}^ 



3 v.. 6 



:V,.((|e|^-3)eMV^,{I±-P±H 
= ±ife • £; ± {I± - P±}«) • E 

+ (i(|^|2-3)Ml/2,L±^i + ^±(^.,u)). 



Furthermore, the system of the high-order moments ([30|) - ([32|) becomes 



(38) dt[@u{{l± ~ P±}") + 2c] + 2d A = e,,(/± + .9±), 
^tQ^A{^± - P±}m) + + d^h, -f V, • {I± - P±}w) 

(39) =e,,(?±-fg±), z^j, 

(40) 9tA,({I± - V±}u) + a,c = + g±). 



The difference between these and dSHl)-® is in ^ since {m^/'^,g±) = 0. 

We conclude this section with some remarks. The derivation of the system pO)) - 
(p2l) or dSgi-dlOl) was initiated by [21 [21], developed for VMB in 37 , and refined in 
[10] by firstly introducing the high-order moment function Q and A. These systems 
play an essential role in the Fourier analysis of the linearized system with general 
microscopic sources as in the case of the Vlasov-Poisson-Boltzmann system [12] 
and even in the study of the hypocoercivity of some linear degenerately dissipative 
kinetic equations [11]. In fact, they are also inspired by the earlier investigation of 
solution spaces without any time derivatives for the well-posedness of the Cauchy 
problem on the pure Boltzmann equation [9] and the Vlasov-Poisson-Boltzmann 
system [13]. 



3. The linearized system with micro sources 

In this section, we are concerned with time-decay properties of solutions to the 
Cauchy problem on the linearized Vlasov-Maxwell-Boltzmann system with micro- 
scopic sources. Specifically, we state the main result in the first subsection, and 
derive a Lyapunov-type inequality for pointwise time-frequency variables in the 
second subsection. Temporal decay rates of the solution and its derivatives in L^- 
norms are obtained under some regularity and integrability conditions on initial 
data and the source terms in the last subsection. 
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3.1. Time-decay properties of solutions. Consider the Cauchy problem on the 

linearized system with a microscopic source g = g(t, x, ^) — [g+, g-]: 

' dtu + i- V^u - E ■ ^M^/2gi =Lu + g, 



(41) 



dtB + V^xE^ 0, 

V, •£;=(Mi/2,w+-M_>, V, -5 = 0, 
I [u,E,B]\t=o = [uo,Eo,Bo], 



where g — {1 — P}g and [uq, Eq, Bq] satisfies the compatibility condition 

(42) V, ■Eo= f M'^\uo,+ - uo.-)dC, V, • = 0. 
For simplicity, we write 

U^[u,E,B], Uo^[ua,Eo,Ba]. 
Formally, the solution to the Cauchy problem (|4ip is denoted by 

(43) U{t) = U\t) + U^\t), 

(44) U'it)=A{t)Uo, U' = [u',E',B'], 

(45) U^'it)= [ A{t- s)[g{s),0,0]ds, f/" = [u", S"], 

where A{t) is the linear solution operator for the Cauchy problem on the linearized 
homogeneous system corresponding to (|¥T|l with g = 0. Notice that U^^ [t) is well- 
defined because [(7(5), 0,0] for any < s < t satisfies the compatibility condition 
(|42| due to the fact that P(?(s) = and hence 



/ m^'^[g+{s)-g^{s)]d^ = 0. 

For brevity, we introduce the norms || • || • with m > and r > 1 given by 

(46) ||C/||2,,„ = + \\[E,B]\\]j^, \\U\\z^ = \\u\\z^ + \\[E,B]\\Lr^, 

for U = [u, E, B], and we set £^ — THP as usual. The main result of this section is 
stated as follows. 

Theorem 3.1. Let 1 < r < 2, I > Q, and let m > be an integer. Suppose that 
(|42]) and Pg = hold. Assume that U is defined in (gSj, (g!]) and (|45|) as the 
solution to the Cauchy problem (I4ip . Then, the first part corresponding to the 
solution of the linearized homogeneous system satisfies 

(47) \\Wl^U'{t)\\c2 < C{1 + t)-|(^-i)-t \\Uo\U + C{1 + t)-^\W^:+'Uo\\c^, 

for any t > 0, and the second part U^^ corresponding to the solution of the linearized 
nonhomogeneous system with vanishing initial data satisfies 



(48) \\^TU"ml^ <C I {1+t- s)-3(^-^)-™||^-i/2g(,)|||^d, 

+ C f\l + t - s)-'\\i^-^^^W'^+^g(s)fds, 



for any t > 0. Here, if £ is not integer, Vf. is regarded as the fractional spatial 
derivative in terms of the Fourier transform. 
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The above theorem shows that solutions to the hnearized homogeneous Vlasov- 
Maxwell-Bohzmann system decays in time with exphcit rates whenever the initial 
data has enough integrability and regularity. In fact, in (1471) . the first term on the 
r.h.s. is generated by the lower-frequency part of solutions for which time rates are 
consistent with those of the classical heat equation, while the second term results 
from the high-frequency part of solutions for which higher regularity of initial data 
implies faster time decay rates. The Maxwell equations of the electromagnetic field 
are essentially responsible for temporal rates with extra initial regularity. This 
kind of phenomenon does not happen to the case of the Vlasov-Poisson-Boltzmann 
system [12], where the second term on the r.h.s. of (|47)) disappears. On the other 
hand, this behavior is consistent with the Vlasov-Maxwell-Boltzmann system on the 
torus 4r. In fact, from the proof of Theorem 13 . II later . it is easy to see that the 
first term on the r.h.s. of (|47|) should disappear if the spatial domain is the torus, 
because this term results from the estimates on the low frequency part of solutions 
and it is also noticed that the frequency variable k takes the discrete values for the 
torus case. Thus, our proof for Theorem 13. II can provide a better understanding of 
the time-decay property of the Vlasov-Maxwell-Boltzmann system at the linearized 
level for both the whole space M.^ and the torus cases. 

The strategy for the proof of Theorem 13. II is to construct a time- frequency Lya- 
punov functional £{t, k) corresponding to the Fourier transform of the system (j4ip 
such that the functional is not only equivalent with the energy-type norm 

\\uit,k)\\l, + \Eit,k)\^ + \B{t,k)\^, 

but also its dissipation rate can be characterized by the functional itself; see Theo- 
rem l3.5l Indeed, the dissipation rate is proportional to the above energy-type norm 
with coefficient 

(l+|fc|2)2- 

Once the Lyapunov-type inequality oi£{t, k) is obtained for pointwise time-frequency 
variables, time-decay of solutions in the physical space follows from the analysis of 
the frequency integration over low and high frequency domains for which p{k) has 
the different pointwise behaviors. 

3.2. A time-frequency Lyapunov inequality. In this subsection, we shall con- 
struct the desired time-frequency Lyapunov functional as motioned before. The 
proof will be carried out along the similar line as in |12| , but additional efforts need 
to be made to take care of the weak dissipation of the electromagnetic field. 

3.2.1. Estimate on the micro dissipation. The first step for the construction of the 
time-frequency Lyapunov functional is to obtain the micro dissipation on the basis 
of the coercivity property (jlip of — L. For simplicity, here and in the sequel, write 
G= (CMi/^u+-w_). Then, 

G - ([e, -^]Mi/2, {I - V}u) = {I - P}u • 

Thus, (|ll])i-(|ll])4 also reads 

' dtU + S,- V^u - E ■ ^M^/^qi =Lu + g, 
dtE-V,xB = -G, 
dtB + V:r X E = 0, 
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Taking the Fourier transform in x gives 

' dtu + i£,- ku- E ■ ^M'^/'^qi =Lu + g, 



(49) 



dtE~ikx B = -G, 
dtB + ikx 0, 
^ ik ■ E = — a^, k ■ B — 0. 



\E\ 



\B\ 



Then equation p9]) i imphes 

iStllull^^ - Re / (Lu I - Re(^ | G) = Re / (g | li)^^ 

We now use the vector identity 

{~ik X B \ E) + {ik X E \ B) ^ 2iRe {k x E \ B), 
to observe from (|49|)2-(l49l)3 that 

^dm^ + \B\')+MG\E)=0. 

Since {E \ G) and (G | E) have the same real part, taking a summation of these 
two equahties gives 

- Re / (Lu I = Re / {g \ u)dt 

JR3 JR3 

From (fTT]) and g = {I — P}.g, one has 

(50) d, (\\u\\l. + \E\' + \B\') + A / - P}u\'d( < C\W-'/'g\\l. 

^ « ^ JR3 « 

where we have used Cauchy's inequahty in the form 

Re/ {gim^Ref {g \ {I - P}u)d^ < %-'/^9\\l, + X\W'/^I- P}i 

for a properly small constant < A < Aq. 

Here, we remark that equation (1501) is the main estimate for the construction 
of the time- frequency Lyapunov functional £{t,k). However, notice that for this 
time, the macroscopic part Pu and i?, B are not included in the dissipation rate of 
([50|) . Next, based on the macroscopic balance laws and high-order moment equation 
obtained in Section [21 we shall introduce some interactive functional to capture the 
rest of the dissipation rate related to Pu and i?, B. 

3.2.2. Estimate on the macro dissipation. Let us apply those computations in Sec- 
tion [2] to the system (|4T|) i-(|4T]1a. Taking the mean value of every two equations 
with ± sign for ([26l), dM]) and noticing Vg = 0, one has 



(51) 



dt 



V, -6 = 0, 



\ 1 3 

2c)+-Y,d,Q,,{{l-V}u-[lA]) = Q, 



dth + d, 

atc + iv, • 6 + ^ ^ 9,A,({I - V}u ■ [1, 1]) = 0, 
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for 1 < I < 3, where moment functions &{■) and A(-) are defined in ([29]) . and we 
used the foUowing facts 

(Mi/2,g^)^(([l,0] + [0,f])Mi/2,g)^0, 

(C,Mi/2,5+ + = (K„e.]Mi/2,5) =0, f < z < 3, 

due to Fg = and hkewise for Lti = [L_|-M,L_ti] due to PLu ~ 0. Similarly, it 
follows from (1301), (EH) and (El that 



(52) 



dt 



ie.,({i-p>.[i,i]) + 2c% 
1 



+a,6j + djh - + + (-9+ + -9-))' 

iatA,({I - P}.. • [1, 1]) + d.c = iA,((Z+ + Z_) + (.g+ + .g_)). 



for 1 < i,j < 3, where Z± is still defined in ([55]) . and (5.^ denotes as usual the 
Kronecker delta. 

Lemma 3.2. There is a time- frequency functional E^^{t,k) defined by 



(53) 



1+ |fc|2 
t2 



ij = l 
3 



l + |fc|2^V^-' 2 



TO^/i fwo properly chosen constants < K2 ^ ki ^ 1 SMc/i i/iai 
A|fcp 



(54) 9tRe£:W(t,fc) 



1 



\b\ 



<C(||{I-P}zi||i. 



holds for any t > and fc G R'^ . 

Proof. The proof can be found in Lemma 4.1]. Although [TTl Lemma 4.1] 
studies the pure Boltzmann equation without force terms our case is directly similar. 
It follows from (|5T|) and ([5^ . which are otherwise not used in the sequel. □ 

In order to further obtain the dissipation rate related to a± from the formula 



\a+ ~ a_| 



2 2' 

we need to consider the dissipation of a+ — a_. For that, taking difference of 
two equations with ± sign for (|26l) . (I?7l) and also noticing Pg = which implies 
{M^/'^,g±) = 0, one has 

(55) dt{a+-a-) + V,-G = 0, 

(56) dtG + V,(a+ - a_) -2E + y,- e({I - P}u ■ gi) 

= (K,-e]Mi/2,g + L{I-P}^.). 
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Note that here and hereafter (V^: • Q)^ (•) = diQij{-). Together with 
(57) V^-£; = a+-a_, 

one has the following lemma. 

Lemma 3.3. For any t > and fc e R'^, it holds that 



(58) + A|a;^_ P < C(||{I - P}^.||i. + 

Proof. In fact, taking the Fourier transform in x for (|55p . ([5^)1 and ([57]) gives 
9t(a+ — a_) + ifc • G = 0, 

^tG + i/c(a^j;'^_) -2E + ik- e({I - P}u • gi) 

= (K,-^]Mi/2,g + L{I-P}zi), 

^ zfc • ii' = a+ — a_. 
On one hand, notice from (15^ that 

(ifc(a![r^_) - I ifc(a^r^_)) (|fcp + 2)|a;:'^_p. 
On the other hand, it follows from ([5^ ? that 

{ik(a^ — a_) — 2i? | ik{a^ — a_)) 

= (-a^G - ik ■ e({I - P}u • -e]Mi/2, 5 + L{I - V}u) \ tk{a^^_)) 

= -dt{G I i/c(a:;^_)) + (G I ikdtia'^^^)) 

-~{ik ■ e({I - P}u ■ gi)-(K, -eiM^/^ 5 + L{I - P}w) I ^^(a;^^-)). 
Combining the above two equations and using (|59Pi . one has 



(59) 



dt{G I zA:(a+-a_)) + (2+|A:n |a+-a_p 

= (G I fcfc-G)-(«fc-e({I-P}M-gi)-([^,-<e]Mi/^5 + L{I-P}M) | ik{a^^^)). 
It follows using Cauchy's inequality that 

5tRe(G I ifc(a^;^_)) + A(l + \k\^)\a:^^-\^ 

< \k ■ G|2 + c\k ■ e({i - p}u ■ qi)f+cm, -m'^',9)f 

+G|(K,-e]Mi/2,L{I-P}7i)|2 

<G(l + |fcni|{I-P}u||i.+G||^-V2^||i.. 

Therefore, ([55)) holds by further dividing the above estimate by 1 + j/cp. Lemma 
is proved. □ 



Notice that from ([57)). the dissipation rate in ((58)) can be rewritten as 

|a;r^_P = \k-E\^ = \k\^\k-E\. 

Here and in the sequel, we always denote k = k/\k\ for |fc| 7^ 0. Therefore, for this 
time, the dissipation rate for A; x ^ and k x B is still not included. Actually, they 
can be recovered from the Maxwell equations of the electromagnetic field as well as 
the evolution equation ()56p of the linear coupling term G. 
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3.2.3. Estimate on the electromagnetic dissipation. As mentioned before, we now 
devote ourselves to obtaining the dissipation rate related to kx E and A: x _B in the 
following 

Lemma 3.4. For any t > and k gM.^ it holds that 
(gQ) f Re(-^fc X ^ I ^) - |fcpRe(G \ E) \ \k x B\^ 



(1 + |A;|2)2 J (1 + |A;|2)2 

Proof. Recall the Fourier transform in x for the Maxwell system 

( dtE-ikx B = -G, 
(61) <^ dtB + ik X E = 0, 

( ik ■ E = a+ — a-, k ■ B = 0. 

It follows that 

\k X Bp = {ikx B\ dfE + G) 

= dt{ik X B\E)-{ikx dtB \ E) + {ik x B \ G), 
where it further holds that 

- {ik X dtB I E) = -{k x{kxE)\E)^\kx E\^. 
Hence, one has the identity 

dt{-ik X B\E) + \kx B\^ ^\kx E\^ + {ik x B \ G), 
which using Cauchy's inequality, implies that 

( Re{-ikxB\E) \ \kxB\' ^ ^ ~ 

To control the first term on the r.h.s. of (15^ . one can again use ([5^ 9 together 
with ^3 and then to get 

- dt{G \ E) + \k ■ E\^ + 2\E\^ 

= -{G I dtE) + {ik ■ e({I - P}u ■ q,y{[t -^]Mi/2, g + L{I- P}u) I E) 
= |Gp -{G\ikxB) + {ik ■ e({I - Pju ■ qi) I ^) 

-(([e,-e]Mi/2,.9 + L{I-PH 1^), 

which implies 

-5tRe{G|£^) + |fc-^P+2|£:|2 

< e{\B\' + \Ef)+^ [{1 + \kf )\Gf + \kfmi P}u ■ q,)f 



-- (m, -e]Mi/2, + G|(K, -e]Mi/2, L{I - P} 



< 6(|Bp + \E\') + ^ [(1 + - Pj^lli. + 
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for a constant < e < 1 to be chosen later. Notice = \k x B\'^ due to 

k ■ B = 0. Then, further muhiplymg the above inequahty by|fcp/(l + |/cp)^ gives 

|fcpRe(G I E) , I (2-^)|fc|V |2 



\k X C |fcp 

'(i + |fc|2)2 +7TT[^2 

Therefore, (pO)) fohows from taking summation of (p^ . and then choosing 
< e < 1 smaU enough. Then Lemma is proved. □ 



Here, we remark that ahhough the pure homogeneous Maxwell system usu- 
ally preserves the energy, the electromagnetic field [E, B] in the Vlasov-Maxwell- 
Boltzmann system indeed has some kind of weak dissipation which results essen- 
tially from the coupling of [E, B] with the microscopic moment function G. 

3.2.4. Derivation of the time-frequency Lyapunov inequality. Now, we are in a po- 
sition to prove 

Theorem 3.5. Let U = [u, E, B] be the solution to the Cauchy problem (|411) with 
= 0. Then there is a time-frequency functional £{t, k) such that 

(64) E{t,k)^\\u\\l, + \[E,B]\\ 



(65) d,£{t,k) + ^-A^£{t,k) < C\\v-^l^g\,% 



where \ {E, B\f = \E\'^ + \B\'^ and for any t > and keR^ we have 

X\k[ 

Proof. Let 

(66) £{t,k) ll^illi. + + «3Re(4V(t,fc) -^4^)(<,fc)), 

for a constant > to be determined later, where S^^^{t, k) is given by (l53j) and 

("2) 

^int fc) is denoted by 

^(2) u ^^ ^-^^ (G I tk{a^^)) {-tk x B \ E) ^ \k\^G \ E) 
[bt) ^int^t,!^)- (i + |A:|2) + (l + |/fcP)2 

One can fix K3 > small enough such that ([M)) holds true. The rest is to check 
(|S5|). In fact, the linear combination of (ISHl) . ([S3]), ([55)) and ([50]) according to the 
definition (j66p implies 

-^1^1^ (\{p\2 , iT. „ i}\2\ / /^ll,.-l/2-||2 



that is 



|2 



(t, k) + Al|^.l/2{I _ p},^j|2^ + _AI^||p^||2^ + x\k\^\k ■ E\' 
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since one has 

\a±f + \bf + \cf^\\Pu\\l.. 

Noticing further that = |A: x Bp due to the fact that B is divergence free, 
foUows. □ 

3.3. Proof of time-decay of linear solutions. Our proof of Theoreni l3.1l is based 
on Theorem 13.51 and the further analysis of (|65p over the low and high frequency 
domains as follows. An alternative new time-frequency splitting method will be 
discussed at the end. 



Proof of Theorem \S.1\ Define the notations 

{i + m^ 

We rewrite (|55|) . for any t >0 and fc e M"^, as 

dt£{t,k) +p{k)£{t,k) < £g{t,k). 
This implies from the Gronwall inequality, for any t > and fc e E'^, that 



£{t,k)<e~P^''^*£{0,k)+ / e-PW(*-^)£g(s,fc)ds. 

Jo 

Due to notice that for any fixed to > we have 

(69) \\VTu{t)f + \\V7[m,B{t)]f- I \kf^£it,k)dk. 



Now, to first prove (H71) . one can apply with g ~ and hence £g{t,k) = to 
bound £{t, k) as follows 

\k\^"'£{t,k)dk < If + [ ] |fc|2"e-f('=)*£:(0,fc)dfc. 

W|fc|<l J\k\>l) 



Here, notice that for |fc| < 1, 



and for |fc| > 1, 



p{k) > 



4|/c|2' 

With these two estimates we have the upper bound 



/ \k\'^"'£{t,k)dk< [ |A:|2™e-^l'=l'*f(0,fc)dfc 

Jr3 J\k\<l 

+ [ |fc|2"e"^£(0,fc)dfc. 



J\k\>i 

Here, as in [351 [55]; from the Holder and Hausdorff- Young inequalities, the integra- 
tion over |fc| < 1 is bounded as 

|fc|2™e-^l^l'*£(0,fc)dfc < Cil+t)-^^^~i^-"^ (lliiolli + ||[i?o,i?o]||i.) , 
'|fc|<i ^ 

for 1 < r < 2. See also, for instance, |12j . 
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The other integration over |fc| > 1 is estimated as 
|fcp'"e"^£'(0,fc)dfc < 



|fc|>i 



|fc|2'"+2^£'(0,fc)dA; sup TTTT^e 
|fe|>i \k\ 



\i.\2e. 

k >l 



1 At 



Since 



sup^e-W<C.,,(l+0- 



it follows that 



/ |fc| 
J|fc|>i 



2 m 



e W£(0,fc)dfc < C(l + t)-^(||V;"+^Mof + ||V™+^[So,Bo] 



Collecting the above estimates as well as (|69p gives (|T7l). 

Similarly, to prove pS)) . one can apply with J/q = and hence £{0, k) = 
to bound £{t, k) as in the following 



/ \k\^"'E{t,k)dk < [ \k\ 



2m 



dk 



ds. 



Recall Sg{s, k) — C\\v ^^^g{s, fc)||^2- Then, similarly as before, it follows that 



11. 



|^|2mg-p(fc)(t-s)^^^^^ fc)dA: < C(l + t - s)-3(^-^)-™|liy-i/2g(s)|ll 

+ C(l+t-s)-^||^-V2v™+Vs)f. 
The above two estimates together with (|69)) give (|48|) . □ 



Remark 3.6. Our initial proof of the time-decay for linear solutions involved a new 
time-frequency splitting method. The key point was to consider the sets 

{t>A|A:|}, {t<X\k\}, 

separately in the energy estimates. In this approach, when t is smaller than \k\ one 
can directly pay for the time decay by placing additional regularity on the initial 
data. And when t is larger than there is a gain to be exploited using the time 
weighted estimates. This approach costs a few extra pages, it could however be 
more robust for other systems in which derivative loss is present. 



4. Energy estimates on the nonlinear system 

The goal of this section is to prove some nonlinear energy estimates, specifically 
(fT5| in Proposition 11.11 and ((T6)) . (|17l) in Theorem 11.21 Combining the time-decay 
property of the linearized system studied in Section |3l these nonlinear energy esti- 
mates can be used to deduce the time rates of the corresponding energy functionals 
in ([T6|) and (fT7|) . which will be shown in the next section. 
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4.1. Dissipation on the electromagnetic field. In this subsection, we prove 
Proposition 11.11 It suffices to prove (fTSj) due to |37l Theorem 1]. In fact, from 
the proof of Theorem 1 in [37], there is £N{t) such that as long as £n{0) is small 
enough, one has 

1<|q|<JV 

(70) + \w^^d^{i-p}uit)r+\\\Eit)r<o, 

\a\ + \l3\<N 

for any < > 0. Recall the definition (fT4|) of ^^^{t) = IJjv.oCi), the rest is to prove 
that the additional term 

l<\a\<N-l 

can be included in the dissipation rate in (1701) . Due to (37. Lemma 6] and the 
uniform-in-time smallness of £jv(^) for solutions, 

(71) Y \\d''Er<c\\{i-p}ur+c Y 

|Q|<Ar-l l<|a|<A' 

This estimate can also be obtained from the following macroscopic balance law 
deduced from the difference of ([55]) with ± sign: 

dM, -e]Mi/2, {I - P}u) + V,(a+ - a_) - 2E + W, ■ e({I - P}u ■ q,) 

= E{a+ +a-) + 2bx B+ ([C, -^JM^/^ Lu + r{u, u)), 

that is, 

2E = dt{[^, -e]Mi/2, {I - P}u) + V,(a+ - a_) + V, • e({I - P}u ■ q^) 

- E{a+ + a^) - 2b X B ~ ([^, -^JM^/^ L{I - P}u + T{u, u)). 

Since E is the sum of some zero-order microscopic term of {I — P}u, first-order 
derivatives and quadratically nonlinear terms, (1711) follows from the above repre- 
sentation of E, and it will be also refined in Lemma [4.21 later on. The dissipation 
estimate on B{t) follows from that of E{t) and the Maxwell system (I7])-(I5]). Recall 

.^2. dtE-V,xB = -([C, -e]Mi/2, {I - P}^), 

^ ' dtB + ^^x E ^Q. 

Take a — (ao, ai, a2, as) with 1 < |a| < A^ — 1. If ao > 0, one can use 
(73) <9"B = d"'dtB = -d"'v^ x E, 

for some a' with \a'\ ~ \a\ — 1 to obtain 

\\d"Bf < Y \\d°'Ef. 

l<\a\<N-l 

Otherwise, when aa — 0, 

d°'B = = -9"'a,A-i(V:r x (v^ x b)) 

= -9"'a.A-iV, X [dtE + ([^, -e]Ml/^ {I - P}u)i 
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for some ct and 1 < « < 3 with |a'| ~ \o\ — 1, where we used the vector identity 
Aaji? ~ ~V X X (Vx X B) since -B is divergence free. Notice that SiA^T^Vx is 
bounded from to itself for any 1 < p < oo. Thus, in the case when ao = 0, one 
has 

||9"S||2<C J2 \\d''Ef + C J2 ll^"{I-P}"ll'- 

l<|Q|<Ar-l |q|<7V-2 

Therefore, taking summation of estimates on B over 1 < \a\ < N — 1 gives 
(74) l|9"S||2<C \\d"Ef + C Y 

l<|a|<Af-l 1<|q|<JV-1 |Q|<Ar-2 

Then p5|) follows from the proper linear combination of ([Tn]) together with ([7T|) 
and (Ell). □ 



4.2. Velocity-weighted energy estimates. In this subsection, we are concerned 
with velocity-weighted energy estimates in the form of (|16p for solutions to the 
nonlinear system (jH])-®. Recall the first part of Theorem II. 2 1 in the following: 

Claim: For any m — 0,1,2,..., there are £N,m{t) and I?Ar.m(t) such that if 
^^Af,in-i(0) is sufficiently small and £^ ^(0) is finite then 

(75) i^N,m{i) + A2?A'.m(i) < 0, 

at 

holds for any t > 0, where £n,-i{Q) — £n{Q) is set and A may depend on m. 

The claim above can be proved by induction on m > 0. First, it is obvious that 
due to Proposition 1 there is £Nfi{t) such that (1751) holds for m = if £n,q{0) is 
sufficiently small. Now, supposing that the claim is true for some m > 0, we shall 
prove that it also holds for m + 1. For that, assume that £N.m{0) is sufficiently 
small. Since £N.m~i{0) < C£N,m{0), then f7v.m-i(0) is also sufficiently small. It 
then follows from the induction assumption that 

(76) ^£N,rnit)+XVN,mit) <0. 

dt 

This implies that £N,m{t) is non-increasing in t, and hence £N.m{t) is sufficiently 
small uniformly in time. Now, our goal is to prove that there are £N,m+iit), 
T^N,m+iit) such that if £N,m+iiO) is finite, then 

(77) jv,™+i (t) + AI?jv,m+i (i) < 0, 

for any t>0. We directly define I'Af.m+iCO by p^ . Then, it remains to construct 
^N.m+iit)- In fact, we carry it out along the lines of HH Subsection 4.1] and [57] . 

Step 1. We split the solution u to equation ([6]) into u = Pu + {I — P}u and take 
{I — P} of the resulting equation to obtain 

(78) dt{l -V}u + i- V,{I - V}u + q{E + ^xB)- - P}u 

= L{I - P}u + r{u, u) + ^E- e{I - P}u + {I - P}iE ■ CM'/^qi) 

- {I - P}($ • V^Pu + q{E + ix B)- V^Pu - ■ ^Pu) 
+ P(e • V,{I - P}u + q{E + ^xB)- VJI - P}u - ^E ■ e{I - P}u). 
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Multiplying the above equation by — P}m and integrating in x,^: 

(79) - P}"ll' - K'+'L{I - P}n, {I - P}u) - Ti + T.+T^, 

where Ti = (i/'"+ir(u, m), {I - P}u), and 

r2 = (1^; • e{I - P}u + P}{E ■ eMi/2gi), - P}^, 

{I - P}(C • VxPu + q{E + ^xB)- V^Pu - ■ ^Pu), - P}u 

+ (P(e • V4I - P}u + q{E + (xB)- - P}u), - P}u) 

^(p(^E-ai-P}u),iy'''+'{I-P}u 



and 

Ts = (-(7(i? + ^ X S) • V^I - P}u), - P}">- 

Here and in the sequel, for simplicity of notations, we also use (•, •) to denote the 
inner product over ^ without any confusion. We will estimate each of the three 
terms in (17^1) . 

These estimates rely on three important observations. The first observation is 
that from |37i Lemma 7] (except for the presence of the momentum weight) and 
also [16| Lemma 3.3] (to handle the momentum weight) we have 

\ri\<c\£ {t)VN,m+l{t). 

The second key observation is that 

-(i/^+iLjl - P}m, {I - P}u) > XWiy'^il - P}uf - Ca||{I - P}uf. 

This follows from L = —h' + K, and the standard compact interpolation estimate 
for K: \{v"'+^KfJ) \ < rjWi^"^ f\\^ + C,j\\f\\^ , which holds for any small tj > 0. 

The third useful observation is the Sobolcv embedding trick used for instance in 
[57] . Specifically we combine the L^(M'^) Sobolev inequality for gradients with the 
embedding W^i'6(]g3) ^ l°°(R3) to obtain 

(80) esssup^gK3 / \u{x,^)\^d^ <C \\u\\l^r^,.d^ 

<C [ ||V,u||^WR,,)de = C||V,u||^i. 

This allows us to control all of the cubic terms in r2 without derivatives. We will 
also use this estimate when only the spatial variables M.^ are present. 
Using Cauchy-Schwarz, we easily obtain the bound 

|r2| < c||ii;|Uoo(R3)||i.^{i - p}u\\' + c\\E\\\\,.Hi - p}u\\ 

+ C\\\/xPu\\\\,.Hl-P}u\\+C\\[E,B]\\L^^^sJ\Pu\\\\ui{I-P}u\\ 
+ C\\iy-2Wx{I-P}u\\\\iy-^{I-P}u\\ 
+ C\\[E,B]\\Lo.^tLl)\W'^^di-P}u\\\W'^{I-P}u\\ 

+ \\EU^^^sJ\,.i{I~P}u\\\ 
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In several places in the estimate above, the velocity growth, ^, is absorbed by the 
projections P. Then it follows from this last estimate and (j80| that 



\T2\ < CVN,mit) + Cyj£N,m{t)VN,„r+l{t). 

To estimate T^, notice 

V^{I - V}u = {'^^v-'^)v'^{l - V]u + v-'^V^iv'^il ~ P}u), 
where the second term on the right contributes nothing into Vj,. Then, one has 

Ta = {-q{E + e X B) • {V ^v-'^)v'^ {1 - P}u, i/^+^jl - P}u) 
.q{m + l){E + £,x B) 



— - V}u,v—{1 - P} 



u) 



' 2v 

Note that we have used v{C) (1 + and |Vjj^(0| < C. We collect the 

estimates in this section to achieve the final estimate of 

(81) i^|lz.^{I-PHp + A||^^{I-PMP 

< CI?Ar,m(i) + £N,m{t)'DN,m+l{t). 

This completes the first step in our proof of (|75p . 

Step 2. Let 1 < |q;| < A^. Applying d" to ([5]) and writing L = —v + K, one has 

dtd"u + ^ • Va;9"u + + C X S) • V^a^u + vd'^u 

= Kd"u + a"r(u, u) + ^d"'{E ■ ^u) + d"E ■ ^M^^^qi 

-qld",{E + ^xB)-V^ju, 

where |-, •] denotes the usual commutator. Multiplying the above equation by 
and taking integrations in x,£, one has 

^l\\^^d-uf + IW'^d-uf = .h + .h. 
In this expression we have used 



Now using the estimates from Step 1, Cauchy-Schwarz, and the Sobolev embedding 
([SO)) , since \a\ > 1, one has the estimates 

I Jll + 1^21 < CVN,m{t) + C^j£N.m{t)VN,m+l{t). 

We collect these estimates, and take summation over 1 < |a| < N, to obtain 
ld_ 
2di 



(82) ~ E \W^d-uf + X Y: IW'^d-^l 



1<|q|<JV l<\a\<N 



This is the main energy inequality in the second step in our proof of ([75 
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Step 3. Let |a| + |/3| < N with > 1. Applying to ^ and writing L = -v+K, 
we observe that u — — P}w satisfies 

(83) dtV + S,- V^v + q{E + i-B)- V^v + lyy ^ h + h + h, 

where li, i — 1,2,3, are defined by 
/i = dpK {d°'{I - P}u) + c'|r(u, u) 

• c{i - P}u) + a^i - P}(a"i? • cmi/^^,), 

h = - P}(e • V^Pu + q{E + ^xB)- V^Pu - ■ £Pu) 

+9|P(^ • V,{I - P}u + q{E + ^xB)- V^I - P}n - |i? • ^{I - P}u), 

and 

/3 = -[9/3, e • V,19"{I - P}u - la^, - P}7. 

-qldl{E + ^xB)-y^l{l-P}u. 

Muhiplying (|83|) by i^^^^w, and integrating over x,^ we have 



Above I,- = (Jj, i/^+iw) for j = 1, 2, 3 and Ia = {{E + ^ x B) ■ V^i/^+i, \qv\^) . 
Exactly as in the analogous estimate in Step 2, we have 



h < Cy £N,rn{t)VN^rn+l{t)■ 
FTOm, for instance, [23l Lemma 2.1], we see that < C so that 

For this we used Cauchy-Schwarz and the Sobolev embedding, as usual. 
We may use the estimate such as [131 Lemma 2.2] to see that Vry > 0: 

\{dpK{d°'{l~P}u),v"'+^d'^{l-P}u)\<f^ \\v'^d'l,{l~P}u\ 

\P'\=\P\ 



Cr,\\v- 



P}i 



Strictly speaking, in [23l Lemma 2.2] there is no velocity weight: However 
this can be added to the proof directly without difficulty. Furthermore, 

I {d$T{u, u), - P}u) I < C^£NMt)'^N,m+i{t). 

This is the content of [371 Lemma 7], the same comment applies for the weight. 
With these estimates, similar to the estimates in Step 2, we have that 



h 



h 



l/3'l = l/31 



■-dZ{i-P}z 
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Since 77 > can be taken arbitrarily small, we add together each of these estimates 
and hnthcr take a summation over I a I + 1/31 < N with 1/31 > 1 to obtain 



|Q| + |^i|<Af \a\ + \0\<N 

I/3|>1 I/3|>1 



nil' 



where Ca,p are some positive constants. This is the third and final estimate which 
we need to prove ([75|) . 

Now, let us define 

fAr,m + l(t) -fjV,m(t) + '«l|k'^{I-P}w|P + Kl ^ W^'^d'^uf 

l<\a\<N 

+«2 E c^./jIk'^a^li-P}"!!', 

|a| + |/3|<JV 
I/3|>1 

for properly small constants < K2 ^ ki <C 1 to be chosen later. Notice that for 
< Ki, K2 < 1, (fT2|) holds true for m + 1. By letting < K2 ^ ki <C 1 be small 
enough, the sum of ([75]) and (ISTjl XKi . (|82p xKi . ([M|) xk9 implies that there is a 
small enough constant A > such that 

■^£N,m+l{t) + AI?JV,m+l(i) < fAf,m(i)2?A',m+l(t)- 

Recall that £N.m{i) is sufficiently small uniformly in time and £Ar.m+i(0) is finite. 
Then, it follows that for any t > 0, £N,m+i{t) is finite and satisfies (|77p . Hence, 
the claim is true for all m > 0. We have shown (1751) and (1161). □ 



4.3. High-order energy estimates. In this subsection, we consider the proof of 
the second part of Theorem ll.2l The goal is to construct a high-order instant energy 
functional £^{t) satisfying the energy inequality (IT71) if £n{^) is sufficiently small. 
For that, we suppose that £n{^) is sufficiently small through this subsection. Due 
to (|15p . £N{t) is also sufficiently small uniformly in time. Recall also the definition 
(m of PjvW- 



Step 1. From the system ([H])-®, the usual energy as in [37] gives 

1 d 
2dt 



l<|a|<Ar l<|a|<Ar 

<c(^£N{t) + ^/£^)vN{t). 

Multiply equation ((75)) by {I — P}w, integrate it in x, ^ and then use ^ to addi- 
tionally obtain 

(86) ~ P}u\f+\\Ef) + X\W'/'{I - P}uf 



< 



[ E-S/^X Bdx + C||t/l/2y^py||2 _^ Cyj£N{t)VN{t). 
JK3 
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Furthermore, similar to (j84|) from the energy estimate on ([78|) . one has 

1 d 
2di 



(87) It, E c^Ad^{i-p}ur + x Ik^5^{i-PMP 



|Q| + |^i|<Ar \a\ + \l3\<N 

I/3|>1 \I3\>1 



<c \W^d''{i-p}uf + c Y (ll3"v,P7/|p + ||a"£;|p) 

Q|<Af |q|<JV-1 



where Ca^p are some positive constants. For completeness, we give the proof of the 
above two energy inequalities as follows. 

Proof of dHU) and §7i): We first prove In fact, ([78]) implies 

^^^^ ^1"^^ " + ^"^^^ ^^""^ ^^"^ 

= ({I - P}{E ■ eMi/2^i), {I - P}u) + {-{I - P}(e • V,Pu), {I - P}u) 

+ (r{u, u) + ^E- e{i - P}u, {I - P}m) 

+ (-{I - PjiqiE + e X B) • V^Pu - ■ SPu), {I - P}u), 

where both the third and fourth terms on the right-hand side are bounded by 
C ^j£N{t)'DN{t)^ and the second term on the right is bounded by 

ry||l.l/2{I_p}^||2^^||^l/2y^p„||2^ 

■n 

for any small constant 77 > 0. To estimate the first term on the right-hand side, we 
use equation ([7]) of _E so as to re- write it by 

({I - P}(£; • eM^/^qi), {I - P)u) ^ I E- ([e, -e]Mi/2, {I - P)u)dx 

= / E-{^dtE + y,xB)dx^-~\\E\\'' + [ E-V.,^Bdx. 
Jr3 I at Jjjs 



Therefore, (|86|) follows by plugging the above estimates into (|88l) . We now turn to 
the proof of As in the proof of dH]), take a, /S with |a| + \l3\<N and \l3\ > 1. 
Applying to ([75)1 and integrating it in x,^ gives 

(89) Ij^imi P}nf + \W'/'d^{I - P}uf = ^(/., - P}u), 

i=l 

where li, i = 1,2,3, are the same as those three terms on the right-hand side of 
(|83)) . The right-hand side terms of ((89)) can be estimated as follows: 

{h,d^{I-P}u)<,j Y \W'/'dp{I^P}u\f 

l/9'l = l/3| 



+ - P}u||2 + Wd'^EW^) + C^£N{t)VN{t), 

(/2,a^{i-p}^)<77l|9^{i-PHr 

+ c^(||i^i/2a"v,Pu||2 + |ia"v,{i - p}uf) + c^/s^v^it), 
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and 



{h,dUi-P}u)<vm{i-P}u\\ 



Cr, J2 \W'^^d^:{I-P}uf + C^lMt)'DN{t), 



\a'\ + \p'\<N 
l/3'l<l/3|-l 

for a constant 77 > small enough. Therefore, ((87)) follows from multiplying (l89l) by 
properly chosen large constants Ca^p > and taking summation over {\a\ + |/3| < 
iV, > 1}. This completes the proof of ([Ml) and (gTl). □ 

Step 2. As in the linearized case, for simplicity, we still denote 

G = (K, ~e]Mi/2, {I _ P}u) = {I _ p}u . qi). 

Corresponding to (|5ip and (|52p for the linearized version, one has from the balance 
laws (|35l) - (|37l) and the high-order moment equations (|38l) - (|40l) in the nonlinear case 
that 



(90) 



a+ + a- 



dt 



0, 



2c)+lY.d,Q,,{{I-'P}u-[l,l]) 



2 



3 



and 



-e„({i-p}ii- [i,i]) + 2c5. 



(91) { +dA + d,h - ^e,,((Z+ + Z_) + (.9+ + .g_)), 
^ iatA,({I - P}.. • [1, 1]) + d,c = iA,((Z+ + Z_) + (5+ + .g-)), 

for 1 < i, j < 3, where l± is still defined in (p3)) and g is defined in (p4l) . Further- 
more, by using the third equation of (|M|) to replace dtc in the first equation of (pij) . 
one has 

(92) i9te,,({i - p}7i • [1, 1]) + d,h, + - ^%v, • 

-^^^jV,-A({I-P}7..[l,l]) 

= ^0..((^+ + + (.9+ + .9-))-^%^ • i?. 

Similarly, corresponding to (ISSI) and ([55]) , it follows from the balance laws ([55]) and 
(PSl) that 

(93) 9t(a+-a_) + V:r-G' = 0, 

5tG + V,(a+ - a_) - 2^; + V,. • e({I - V}u ■ qi) 

(94) = E{a+ + a_) +2bx -^JM^/^ Lm + r(w, u)). 
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We also recall 

(95) yx-E = a+-a^. 
Lemma 4.1. One has the following four estimates 

(96) W^Ab,cW<C \\d"{I-P}uf + C£N{t)VN{t), 

|q|<7V-1 |a|<A' 



(o-i- + aJ)dx 



l<|a|<W-l 



2 



A l|a"V,(a+ + a_) 

l<\a\<N-l 

<C Y \\d"{I-P}u\f + C£N{t)VN{t), 

\a\<N 



(98) 4 E / 9"G-9"V,(a+-a 



a|<Af-l 



— a^)dx 

A E (||9"V,(a+-a_)|p + ||a"(a+-a_)f) 

|a|<Af-l 

<C E ||9"{I-P}u||2 + C£:jv(i)25iv(0, 

\a\<N 



and 



(99) E ||a"atK±a-,6,c]||2 <C Y l|5"V,(a++a-)|P 

|q|<JV-1 |Q|<Af-l 

+ C E \\d''{l-P}uf + CEN{t)VN{t), 

\a\<N 

for any t> 0. 

Proof. First consider (1961) . For a with |a| < — 1, the bound of c follows from 



WdA' < C\\d"dtU{l - P}u ■ [1, 1])||^ + C||a"A,((/+ + /_) + {g+ + .g_))|| 
< C\\d''dt{I - P}uf + C||9"V4I - P}uf + C||9"{I - P}uf 

+C||9"A,(.g±)|p 
<C7 E ||9"{I-P}u|P + C£^(i)I?^(t) 

|q|<JV 

due to the second equation of (|9T|) . and in the same way, the bound of b follows 
from equation ((5^ by noticing the identity 



/ \d"id,b, + djb, - 9%V, • b)\''dx = 2|ia"V,6f + ^|19"V, • bf 

7r3 



— ( 
3 
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Next, (j97| follows from the first two equations of (|90l) . In fact, for a with 1 < |a| < 
— 1, the second equation of (|90l) implies 

i||a"9,(a+ + a_)f 

d°'d,ia+ + a^)d"b,dx + 2 [ d°'d^dt "'"'t"'~ d'^hdx 

+ / a"a,(a+ + a_)a"|-2a,c-i Va,e.,({i-p}u- [1,1]) 

+ +[GxB],}dx, 
and further using the first equation of (|90l) and Cauchy-Schwarz gives 

V / 9"6-a"V,(a+ + a_)da; + A V ||5"V,(a+ + a_)|p 

l<\a\<N~l'"^ l<\a\<N-l 

<C \\d"V,[b,c]f + C J2 \\d''yx{I-P}u\\^ + C£N{t)VN{t) 

l<|a|<7V-l l<|a|<Ar-l 

<C J2 \\d"{I~F}u\\^ + C£N{t)VN{t), 

\a\<N 

where was used in the last inequality. Hence, (1^7)) holds. To prove (|M1) . for a 
with |q:| < ^ - 1, dnH) together with dM]) and ^ yield 

||9"V,(a+-a_)||2 + 2||a"(a+-a_)||2 

/ d°'W^{a+-a^)-d°'[W^{a+-a^)-2E]dx 

dt 7j{3 

+ / 9"V^(a+ - • a" [ - • e({I - P}u ■ qi) + E{a+ + a_) 
+2& X B + ([e, -e]Mi/2, L{I - P}u + r(u, u))]dx, 

which implies (j98p after using Cauchy-Schwarz and taking summation over \a\ < 
N — 1. Finally, it is straightforward to verify ([99|) from (|90| and (|93| as well as 
(|96| . The proof of Lemma 14.11 is complete. □ 



Furthermore, ([7T|) and ([71)1 for the upper bounds of E', B can be refined as 
Lemma 4.2. It holds that 

(100) ^ \\d"Ef+ J2 \\c>"Bf 

la|<Af-l l<|a|<Af-l 

<C ^ ||9"V,Pu|l2 + C ^ ||9"{I-P}u||2 + C£:w(t)I?^(i), 

Q|<Ar-l |q|<A' 

for any t > 0. 
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Proof. One again, we use ([94|). that is, 

(101) 2E = dtG + V^{a+ ~ aJ) + • e({I - P}u • qx) 

- E{a+ + a_) - 26 X B - ([C, -C]M^/^ L{I - P}u + T{u, u)). 

Then, the upper bound for E in (jlOOp fohows directly from the above equation, 
and thus the upper bound of B also holds by (TMl) and the estimate on E. □ 

Step 3. Let us define the interactive high-order instant energy functional jj^^ {t) 
as in [13] by 



(102) £'^-Ut)= I d"b-d"V,{a++a.)dx 

l<|a|<A'-l 

+ V / (9"G-a"V^(a+-a_)dx. 

Notice that 

(103) K,,M<c J2 \\9"Puf + c Il5"{i-PM^. 

l<|Q|<Af |a|<Af-l 



In addition, the proper linear combination of (|96| . ((97|) . (|98p and ((99)) in Lemma 
14.11 implies 

±£%,^^{t)+X Y \\0"[a+±a.,b,c]r + X\\a+-a^r 

l<\a\<N 

< C||V,(a+ + a_)||2 + C Y \\d''{I-P}uf + C£N{t)VN{t). 

\a\<N 

Note that to derive the above inequality we can add ||Vj;(a+ + a-)|p to both sides. 
After further plugging in (jlOOp . one has 

p^,^^it)+X Y l|5"K±a_,6,c]f + A||a+-a_f 

l<|a|<Af 

q|<JV-1 l<|Q|<Ar-l 

< C||V:,Pu|P + C \\d°'{I-P}u\\^- + C£N{t)VN{t). 

\a\<N 

Hence, it further holds that 

(104) |f]^,„tW+A Y ii5"P"iP+A Y mE,Br+x\\Er 

l<|a|<A' l<|Q|<Af-l 

<C\\W^Puf + C Y Wil-Pjuf + C£N{t)VN{t). 

\a\<N 

This is the main energy estimate for £^^^^^{t). 
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Step 4- Notice that ([86]) implies 

(105) ~m - pjur+wEf) + p}ur 

+ C {£N{i) + ^JI^)vN{t) 

for a constant 7y > to be chosen small enough. Now, we are ready to construct 
S'^it). In fact, let us define 

^N{t)= E {Wuf+\\d-[E,B]r)+n,m-^}nr+\\Ef) 

l<\a\<N 

\a\ + \fi\<N 
I/3|>1 

for suitable constants < K3 ^ K2 ^ ki ^ 1 to be determined now, where £^ jj^j {t) 
is given by (|102p . Due to (|103p . one can let < K3 ^ ^2 ^ ki <C 1 be small enough 
such that holds with m — Q and hence £^ {t) is indeed a well-defined high-order 
instant energy functional. In addition, by choosing < K3 <C K2 ^ ^ 1 further 
small enough, the sum of (|85|) . (|105p xKi for 77 > small enough, (|104p xK9 and 
([571) X K3 yields 

+ XD^it) < C(||i^i/2v,Pu||2 + ||V, X B\\') 

at 

+ C(^y/£Nit)+£Nit))VN{t), 

which implies the desired estimate (fT7|) since £n {t) is small enough uniformly in all 
t > 0. This completes the proof of the second part of Theorem 11.21 □ 

5. Time decay for the nonlinear system 

This section is devoted to the proof of Theorem 11.31 and Corollarv ll.4l As men- 
tioned at the beginning of Section 21 the time-decay rates of solutions in some 
energy spaces are obtained by combining the corresponding nonlinear energy esti- 
mates and the linearized time-decay property applied to the lower-order terms. In 
addition, we also will use the time- weighted estimates, as in (|107p . and a bootstrap 
argument on the derivatives to handle the regularity-loss phenomenon in the energy 
dissipation rate due to the degenerately dissipative (regularity loss) property from 
the Maxwell system. Decay rates of solutions in the L°°-norm are deduced by using 
some optimal Sobolev embedding theorem, and hence time rates for L''-norm with 
2 < r < 00 will follow from the normal L'^-L°° interpolation inequality. 

5.1. Decay rates of the full instant energy functional. In what follows, 
we fix an integer m > and suppose that eN+2,mvi defined in ()18|) from initial 
data [uq, Eq, Bq] is sufficiently small. The goal is to prove ([T^ . Since £N.m{0) < 
eAf+2,mvi is sufficiently small, the Lyapunov inequality (1T5|) for £N,m{t) and VN^mit) 
holds for any i > 0. Let 1 < ^ < 2. Multiplying ^ hy {I + tf and then taking 
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time integration over [0,t] gives 

(1 + tYSN.mit) + A / (1 + sYVN,m{s)ds < £N^,niO) +i I (1 + s)'"'^ iV,m(s)ds. 



Recall and (|14l) for definitions of £N.m{t) and I?A',m(t). It follows from 
and dUl) that 



£N,mi-t) < Ci\\Fu{t)f + \\B{t)f + VN+l,^{t)). 

Note that we could take I'Ar+i.m-i(i) above. Then, it follows that 

(l+t)^fw,„(t)+A / {l+sYVN,,n{s)ds 

Jo 

< £n,UO) +C£ f (1 + sY-^ {\\Pu{s)\\^ + ds 



+ Ci (1 + sY-'VN+i.mis)ds. 

JQ 

We use similar estimates from ((TE)) for + 1 and + 2 as above to obtain 

(l+t)^-lfjV+l,mW+A / {l + sY-^VN+l,„^{s)ds 
JQ 

< £N+i.m{0)+C{e~l) f {1 + sY-H\\Fu{s)\\^ + \\B{s)\\^)ds 
Jo 

+ Ci£-1) f (l + s)^-'Pjv+2,™(s)ds, 



and 

£N+2.m{t) + A / VN+2,m{s)ds < f Ar+2.m(0). 

Jo 

It follows by iteration that 

(106) (1 + tY£N,m{t) + A / (1 + sYVN.m{s)ds 

Jo 

< C£N+2,m{0) + Ce f {1 + sY-\\\Fu{s)\\^ + \\B{s)\\^)ds, 
Jo 

due to 1 < ^ < 2. 

On the other hand, define 

(107) £?^„,{t) = sup (1 + s)-^£N,rn{s). 

0<s<t 

From now on let us consider the pointwise time-decay estimate on ||Pu|P + WBW^ 
in terms of f^„(t). Formally, the solution U = [u,E,B] to the Cauchy problem 
(HI) , © , ([S]) , © of the Vlasov-Maxwell-Boltzmann system can be written as 

(108) U{t)=Ait)Uo+ [ A(t - s)[5(s),0,0]ds, 

where A is defined in and g — [g+,g^] is given by p4p : equivalently 
g r(u, u) + ^E-^u-q{E + ^x B)- Vjm. 
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For later use, U (t) can be rewritten as 

Uit)^Ioit)+h{t)+l2{t), 



with 



loit) ^ A{t)Uo, 

hit) = I A(t-s)[5i(s),0,0]ds, 



/2(t) - / A(t-s)[52(s),0,0]ds, 
Jo 

where 

51 = r(u, u) + {l- P} (^^E -Cu-qiE + ^xB)- V^u) , 

52 = P [^E -^u-qiE + ^x B)- V^w) . 

Here, we remark that time- integral terms Ii(t), hit) and the time-integral term in 
Uit) are well-defined since 

(109) / Mi/2[5^^^(s)-5,^_(s)]de = 0, 1 = 1,2, 
for all X e M^. For later use, we will need 

Lemma 5.1. Assume N > 4. For any integer j > 0, there is C > such that 

(110) ||:.-i/VWIIi.(H^)nz, =^^^^v^,iWf,v^(i), 

(111) \\92it)\\L2^Hi)nZ, < CS.^Nit), 

for any t > 0, where j \/ N = max{j, N}. 

Proof. Due to the definition of the quadratically nonlinear function 52 and the fact 
that P is a macroscopic projector, (jllip immediately follows by the definition 
for SNit) = SNfiit) with iV > 4. To prove (ITTO| . 

+ 2\\iy-^/^ (1^; -Cu-qiE + ^x B)-W^ 

+ 2\\u-^/^P (l-B -^u-qiE + ^x B)-V^ 



2 

2 



where it is also straightforward to see that the last two terms on the r.h.s. are 
bounded by C£jyN,iit)£jvNit), while the estimate on the first term follows from 
[TBI Lemma 3.3]. Hence, (IllOp and (|llll) are proved. □ 

We recall the norms from (l46l) for [/ = [u,E, B]. We now apply (|47| with to = 0, 
r = 1 and ^ = 3/2 to /o(i) and /2(i), respectively, to obtain 

\\Ioit)\\c^ < Cil + t)-i\\Uo\\^s^^^, 
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and 

Whmc^ < f \\A{t-s)[g2is),0,0]\\c2ds 
Jo 

ft 

< C / (l+t-s)-3||g2(s)|| 3 ds. 

For l2{t), it further holds from Lemma [01 and the definition (|107p of £^,-^{t) that 

\\h{t)\\c^<c[ {l + t- syi£N{s)ds<C j {l + t- s)-i£N,m{s)ds 
Jo Jo 

_ 3 







Similarly, applying (|48)) with m — 0, r = 1 and £ — 3/2 to Ii{t) due to the fact 
that Pgi — and then using Lemma [STTl one has 



\\Iimh<C {l + t-s)-^\W-'/'g,{s)\\\ 3 ds 
<C {1 + t- syhNAs)£Nis)ds 







<C / (l+t-s) '^£N,„iVlis)£N.niis)ds. 

Jo 

Further using the definition ()107p of again and (fT6|) . it follows that 



||/iWll2:2 <C£:^„(Osup£:w,™vi(r) / (1 + i - s)-^ (i + 
T>0 Jo 

< C{1 + t)-h^^^{t)£N,mMQ)- 

Collecting the estimates on Ii{t) [i — 0, 1, 2) above implies 

2 

(112) wpu{t)r + < c\\u{t)\\i. < wmwh 

i=0 

< C{1 + t)-^\Uo\\l^2nz, + C£N^m^i{0)£^^^{t){l + 

+ C{l + t)-i[£Z^{t)f. 
Now, fix a constant e > close to zero. Taking £ = 3/2 + e in (|106p yields 

(l + t)i+'£jv,™W + A / {l + s)^+'VN,m{s)ds 







< C£N+2,rn{0)+C I \l + s)i+'-' i\\Puis)f + \\B{s)f)ds. 

Jo 

Plugging ()112p into the above inequality gives 

{l+t)i+'£N,m{t) < C£jv+2,m(0) 
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which after dividing by (1 +tY imphes that for any t > 0, 

(1 + t)i£N,m{t) < C (fjV+2,™(0) + ||t/o||^.n2, +^iV,mVl(0)£:^,„(t) + [^K^mW]') • 

We equivalently have that 

f^,™(t) < C (£Ar+2,™(0) + ||C/o||?,.nZ, + £N^.^^m£N.J^ + [^N.M^) ■ 

Since £N+2,miO) + l|C^o||«2n2i < eAr+2,mvi and £N,mvi{0) < eN+2,mvi are suffi- 
ciently small, one has 

,771 : 

which gives the desired time-decay estimate (IT9l) . □ 

5.2. Decay rates of the high-order instant energy functional. In this subsec- 
tion, we shall prove the time-decay estimate (1201) for the high-order instant energy 
functional £%{t). It follows from Theorem 1 1 . 21 that 

(113) ^^£Ut) + APjvW < C(||z.i/2v^Pu(i)||2 + ||V, X i?(i)f ), 
which by definitions ([T^. pi)) of £]i^{t) and I?jv(i), implies 

(114) j^£'^{t) + X£Ut) < C(||i.i/2v^P7.(t)||2 + ||V, X B{t)r) 

+ c J2 \mEit),B{t)]f- 

\a\=N 

Notice that from the system ©-([S]), one can replace all time derivatives of E and 
B by the spatial derivatives of E, B and space-time derivatives of {I — P}u. 
In particular we claim that 

(115) ^ ||9"[i?(t),i?(i)]|P<C II{I-P}5"'^^(0IP 

\a\=N |Q'|=Af-l 

+c j2 mEit),Bit)r. 

\a\=N 

Q0 = 

Proof of claim: This estimate (|115p follows from an induction on the number of time 
derivatives using ([7^ . First suppose \a\ — N and = 1. Then a = (ap, 0, 0, 0)+a' 
with a' = (0, Qfi, Q!2, cks) and |a'| = iV — 1. Now using ([75]) we obtain 

\\d"Bf <C J2 W^^'d^Ef <C J2 

|a| = l l<|a|<Ar 
cto— ao— 

Similarly, using (j72p and also Cauchy-Schwarz, in this initial case we have 

|q| = 1 

q;o=0 

This completes the first inductive step. 

Now we suppose that the lemma is true for |a| = N and aa = m with any 
TO = 1,...,A^— 1. We will show that the lemma is also true for ao = m -f- 1 . Similar 
to the first inductive step, as in ([75)) we have 

d°'B = d°''dtB = -c)"'v^ X E, 



38 



RENJUN DUAN AND R. M. STRAIN 



for some a' with \a'\ = |a| — 1 since ao = to + 1 > 0. Since ag = to, we achieve 

\\d"Bf <C Wd^'d^'Ef. 

\a'\=N-l \a\ = l 
ag=m "0=0 

Now the general estimate for B foUows from the inductive hypothesis. 
For the last term, we use ([72]) i to observe that 



d°'E = d°''dtE = 9"'V:r X S - ([C, -^JM^/^ {I - P}^" u). 
Similarly, using Cauchy-Schwarz we observe that 

a'|=Ar_l |q| = 1 |a'|=Ar-l 
a'a=m "0=0 

This completes the estimate (jll5l) by induction. □ 

Now, taking a suitable linear combination of (|114p and (|113p , we see from (|115p 
that (|114p can be rewritten crudely as 

|f^(t) + A£]^(t)<C7 {Wu{t)f + \\d-[E{t),Bm?)- 

l<\a\<N 
"0=0 

Here we noticed that from for a small k > it holds that 
VN{t)~K Y IIII-Pl^" w(<)f > 0. 

|a'|=JV-l 

Then using the norm (|46p . we have equivalently shown 

(116) + - C'||V,C/(t)||^«-i. 

As before, it follows from Theorem 13.11 and the representation (|108p of U{t) that 

(117) ||V,[/(t)||^„_, <C(l + i)-i||[/o||?,«+3n2, 
+ C / (l + i-s)-t||i/-i/V(s)||'2(^«+3)nZi'^« 



C 



(1 + i - s) 4 ||5r2(s)|li2(^«+3)r^2^ds 



Recall the estimates from Lemma |5. II so that one has 

(118) \\^'^'^9l{ml2(H-^')nZ, <C£N+3,lit)£N+3it), 

(119) ll52(0llL|(Hf+3)nZi < C''^^W+3(0> 

for any t > 0. Since we suppose that ejv+5.1 is sufficiently small, the first part of 
Theorem 11.31 implies 

£N+3,l{t) < CeAr+5,l(l + 0"^- 

Putting (I118p - (|119p together the above estimate into (|117p then gives 

||V,[/(t)||^„-i < Cejv+5.i(l + 0"^- 

Using (|116p . this implies from the Gronwall inequality that (f20|) holds true for any 
t > 0. This completes the proof of Theorem 11.31 □ 
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Remark 5.2. It is easy to see that, using the similar method as for the proof ([T6| . 
(fT7|) without velocity weights in the high-order energy functional £%{t), the Lya- 
punov inequality can be refined as 

for any t > provided that £N.7n-i{0) is sufficiently small, where the integer m > 
is given. Thus, starting from the above inequality as for (|113|) . it can be proved 
that 

holds for <: > if ^^Ar.m-i(O) and eAr+5,1 are sufficiently small. 

5.3. Decay rates of solutions in L^.. In this subsection, we shall prove Corollary 
11.41 It suffices to obtain time rates in (I^TI) and ([^^ . They hold true when r ~ 2 due 
to assumptions of Corollary 11.41 Theorem [O] and definitions p^ . oi £N,m{t) 
and £^j^{t), where (f7T|) is used to estimate 
Then, the rest is to verify 

(120) \\Umz^^\\u{t)\\L.^i^^) + \\[E{t),BmL^ <C{l+t)-i. 
Lemma 5.3. Let Uq — [uq, Eq, Bq] satisfy (|42p . Then, it holds that 

(121) \\A{t)Uo\\z^<C{l + t)-^\\Uo\\z, + \\^lUo\\n-), 
for any t > 0. 

Proof As in dH]), let us set U'it) = A(t)Uo with = [u',E',B']. From the 
Sobolev inequality (43, Proposition 3.8], one has 

(122) \\u'U^^(L^^ < C\\W,u'\\'/^WW\\'/\ 

Notice that the proof of this inequality also follows easily from optimizing ((80)) . For 
the first-order derivative, using (|47p with m = l,r = 1,^ = 5/2, it follows that 

||V,u^(t)|| < C{l+t)-i{\\Uo\\z, + WVlUoWn^). 
For the second-order derivative, similarly letting m = 2,r = 1,^ = 7/2 in (1471) . one 

W'^Wm < c{i + t)-i{\\u,y, + ||vEc/o||«0- 

Putting the above two estimates into (|122l) implies 

\\u'it)\\z^<C{l+t)-li\\Uo\\z, + \\VlUo\\n^). 
Similarly, it holds that 

\\[E\t),B'mz^<C{l+t)-H\\Uo\\z. + \\VlUo\\n^). 
Therefore, follows. □ 

Now, by applying Lemma 15.31 to the representation (|108p of U (t) , one has 

(123) \\U{t)\\z^ < Cil +t)-^\Uo\\n^nz, 

+C / {I + t - s)-^\g{s)\\L2fH^Jnz^ds. 
Jo ' 

Notice from Lemma [5. II that 
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Now since £8,2 in ([T8|) is sufficiently small, it follows from (fTQ]) that 

Putting the above estimates into (|123p leads to (|120p . This completes the proof of 
Corollary O □ 
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